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A recent report that sulfur hydride under pressure is an electron-phonon superconductor with 
a T c of 190 K has been met with much excitement although it is yet to be confirmed. Based on 
several electron-phonon spectral density functions already available from density functional theory, 
we find that the electron-phonon spectrum is near optimum for T c with a particularly large value 
of its characteristic phonon energy cvi n which is due to the small hydrogen mass. We find that the 
thermodynamic universal BCS ratios are near those for Pb and NbaSn. We suggest that optical 
measurements could be a useful tool to establish the existence and nature of the superconductivity in 
this system. Conventional superconductors are in the impurity-dominated dirty limit. By contrast 
sulfur hydride will be in the clean limit because of its large energy gap scale. The AC optical 
conductivity will display distinct and separate signatures of the superconducting gap in the low- 
energy impurity-dominated range of the optical spectrum and additional phonon structures at higher 
energies where the clean limit applies. 

PACS numbers: 74.20.Fg,74.62.Fj,74.25.Gz,74.20.Pq 


The observation^ of a dramatic and sharp drop in resis¬ 
tivity around a temperature of T < 190AT in a sulfur hy¬ 
dride sample under ~ 200 GPa, has led to a flurry of new 
theoretical works^F— . The verification that it is definitely 
superconducting, for example through the observation of 
a Meissner effect, still awaits further experiments, how¬ 
ever, the authors did observe that with the application 
of a magnetic field of 7 Tesla there was a shift, to lower 
temperature, of the onset of the resistivity drop by ~ 10 
K. Moreover, the conductivity below 190 K was found to 
be better than that of the purest metal at low temper¬ 
ature and was two orders of magnitude larger than that 
of pure copper. 

A number of papers^r— based on density functional the¬ 
ory for sulfur hydride have already provided detailed in¬ 
formation on the electron-phonon spectral density func¬ 
tion a 2 F(oj) and have identified H 3 S as the most likely 
phase involved at that pressure. While some quantitative 
differences exist between the various spectral densities 
available, all calculations support the idea of phonon- 
mediated superconductivity in H 3 S with a T c of order 
200 K at ~ 200 GPa. Still it remains essential to 
find ways to directly verify such a possibility. Lack¬ 
ing definitive confirmation, it is instructive to place H 3 S 
within the context of the family of known conventional 
superconductors^. In particular an emphasis on dif¬ 
ferences and commonalities between H 3 S and the well- 
established cases of electron-phonon-mediated supercon¬ 
ductivity can provide valuable insight. It is also impor¬ 
tant to calculate the expected properties of the supercon¬ 
ducting state of H 3 S with the specific aim of providing 
information that can be used by experimentalists to es¬ 
tablish that H 3 S is indeed a superconductor and to dis¬ 
cover the pairing mechanism involved, i.e., the glue that 
binds the Cooper pairs. 

In this Rapid Communications, we demonstrate that 


H 3 S is a very highly optimized electron-phonon super¬ 
conductor for various reasons and we provide a calcu¬ 
lation of the optical conductivity as a test for identi¬ 
fying the superconducting gap and the Holstein struc¬ 
ture for the phonon spectrum. For this purpose, we 
work with the u 2 F(oj) spectrum provided in the liter¬ 
ature. For illustrative purposes we have used the one 
developed by Errea et alA (shown in Fig. IS1I) which in¬ 
cludes anharmonic phonon effects. However, in the Sup¬ 
plemental Material^ (Fig. SI), we also show the spectra 
for the harmonic version from Errea et al. along with 
that from Flores-Livas et alA for comparison [the one 
from Duan et al 2 is also similar but not shown]. Us¬ 
ing the a 2 F(ui), one can calculate the electron-phonon 
mass renormalization parameter, A = f Q [a 2 F(co)/uj\du>, 
the area under the spectrum A = J 0 °° a 2 F(uj)duj and the 
Alien-Dynes characteristic phonon frequenc y 11 : 12 uq n = 
exp[(2/A) / 0 ° o [lii(w)a 2 F(w)/w]da;], which would equal oje 
if the a 2 F(oj) = A5(u> — oje)- For our digitization of 
the published spectrum shown in Fig. ED A = 1.67, 
A = 118.5 meV, and wi n = 122 meV. Comparison of 
these characteristic properties with the other spectra are 
given in Table SI of the Supplemental Material^. [Note, 
these spectra were digitized from the original papers and 
placed on a 2 meV grid. This gives characteristics that 
are in qualitative agreement with the original works if not 
always in perfect quantitative agreement. Such details do 
not impact our results and conclusions presented here. It 
is clear that all of the spectra are qualitatively similar 
Calculating T c with these spectra using the standard s- 
wave electron-phonon Eliashberg equations^, we find that 
for T c = 190 K, the Coulomb pseudopotential pA would 
have to be 0.18 for the anharmonic spectrum and about 
0.38-0.4 for the harmonic spectrum (or conversely, for a 
pA ~ 0.16, the harmonic spectra would give a T c of or¬ 
der 250 K, which has been noted by the authors of these 
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spectra). While a large /z* is unusual, such a value has 
been found for the alkali-doped fullerenesi^ and can be 
taken as an argument for strong electron correlations to 
be important in addition to the electron-phonon interac¬ 
tion. It should be noted that the spectrum of Flores-Livas 
is one that has been developed and used by accounting 
for a nonconstant energy-dependent electronic density of 
states and yet it essentially agrees with the Errea har¬ 
monic spectrum. These details aside and working with 
the anharmonic spectrum for illustration, we argue that 
the a 2 F(u}) for H 3 S is one of the most highly optimized 
among conventional superconductors. 

One test of this concept is to evaluate the functional 
derivative of T c [Ref-^] with respect to this spectrum 
as is shown in Fig. ISTT al (red curve). This derivative is 
peaked at u> ~ 7 ksT c and illustrates that if one could 
move all the spectral weight in the a 2 F(uj) to a delta 
function placed at this frequency, then one would obtain 
the highest T c possible from this spectrum. However, 
the spectrum is distributed around this point and so the 
question remains as to how optimized is it? Using the 
characteristic phonon energy uj\ n as a measure, one finds 
that it is very close to the value of 7T C as shown in the 
figure. Consequently, the spectrum is quite highly opti¬ 
mized already. This is understood from the fact that the 
functional derivative varies quite slowly about its peak 
and much of the weight of the a 2 F(uj) falls in the re¬ 
gion of this slow variation. The shape of the functional 
derivative speaks to the physics that for phonon frequen¬ 
cies that are too low, the phonons do not respond fast 
enough to pair the electrons, but rather act more like 
static impurities. Alternatively, at very high frequencies, 
the ions would wiggle too fast to have any net displace¬ 
ment or polarization to glue the Cooper pairs. As a result 
the optimal frequency is found at some intermediate fre¬ 
quency. A simple physical argument for the maximum to 
occur at w ~ 7 ksT c can be found in Refs^^ 

The question as to how optimal this spectrum is can be 
quantified by examining Fig. ISlI bl where the dimension¬ 
less ratio of T c /A is plotted versus /i*. The line with the 
open circles is the maximum value that T c could take for 
fixed A as a function of /z* and is calculated by placing a 
delta function of weight A at the frequency of the peak in 
its own functional derivative and finding the T c . This is 
a classic curve taken from Leavens^ and the solid black 
dots represent the data for several conventional super¬ 
conductors taken from a tabulation (Table IV) in Ref£. 
The gold star shows where the spectrum in Fig. ISTT a! sits 
and it is just below the optimal line. Indeed, it falls only 
8 % below the theoretical absolute upper bound allowed 
to any spectrum with that area A and /z*. Furthermore, 
the blue curve, shows how T c would vary for this spec¬ 
trum as a function of /z* and it is clear that it would 
always remain close to the upper bound. For /z* = 0, the 
T c = 298 K. 

One may now try to understand from this why the 
T c is so high. It is the unusually large value of A in H 3 S 
which accounts for the large T c found in this system. This 




FIG. 1. (Color online) (a) The anharmonic a 2 F(u>) spectrum 
for H 3 S at 200 GPa given in Ref.—. The functional derivative 
STc/Sa 2 F(lo) vs w for this same spectrum is also shown scaled 
up by a factor of 12. It peaks at uj = 7It bT c (dashed line) 
which is very close to the characteristic phonon frequency 
win (dash-dotted line) for the a 2 F{uj) demonstrating that the 
spectrum is highly optimized for T c . (b) T c /A vs /z* where A is 
the area under the a 2 F(uj) spectrum. The black curve with 
open circles is the theoretical upper limit for T c with fixed 
A and the solid circles are data for several superconductors 
constructed from Table IV of Reh 9 . The gold star is for the 
a 2 F{uj) spectrum in (a) with T c = 190 K and /z* = 0.18. The 
blue curve represents how T c /A would change for varying /z*, 
demonstrating that this spectrum is near optimal. 


sizable A comes about because the electron-ion coupling 
remains large to very high phonon energies. In turn, the 
unusually large value of w>i n is traced to the small hy¬ 
drogen mass and the stiffness of the H-S bond stretching 
mode. Moreover, the lack of an electronic core struc¬ 
ture to the hydrogen (see particularly Ref4) means that 
there is no orthogonalization of the conduction electron 
wavefunctions which in ions with core electrons effec¬ 
tively changes the electron-ion interaction potential into 
a weakened pseudopotential. These two facts combined 
provide the mechanism for a large area A under the spec¬ 
tral density a 2 F(u>) which extends to high energies with 
large magnitude. Further insight into the importance of 
a large wi n can be obtained through an analysis based 
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on the asymptotic formula- for T c of Eliashberg theory 
derived on the assumption of a very large value of A, the 
mass enhancement factor. In this limit T c = 0.182\/Au;i n . 
We see that increasing A by a factor of 2 increases T c by 
~ 40% while increasing wi n by 2 doubles the T c . While 
H 3 S is not in the asymptotic limit of large A, this analysis 
nevertheless provides a useful qualitative picture. While 
A for H 3 S is not particularly large as compared with other 
electron-phonon superconductors, wi n is quite significant 
at about 120 meV. Through the numerical solution of 
the Eliashberg equations for many model spectra, Leav¬ 
ens and Carbotte!! observed that, for A and ji* in the 
range of 1.2 < A < 2.4 and 0.1 < fi* < 0.15, a re¬ 
markably simple relationship between T c and A applied, 
namely T c = 0.148A This formula can be used to obtain 
further understanding of the role various regions of the 
a 2 F(uj ) spectrum play in T c . As an example, the entire 
region below ~ 165 meV in the spectrum of Fig. IS If al 
has an area of 39 meV while the area above is 80 meV, 
contributing, respectively, 67 and 137 K to the T c value. 
It is clear that the hydrogen bending vibrations and most 
importantly the stretching modes are critical for the high 
T c of H 3 S metallic phase. 

In Fig. [S2l we compare the gap ratio calculated for 
sulfur hydride with other conventional superconductors 
on a classic plol— for the gap ratio versus T c /oj\ n . We 
have done this calculation for both the anharmonic spec¬ 
trum (gold star) and the two harmonic spectra shown in 
the Supplemental Material^ (with T c = 190 K). Clearly, 
sulfur hydride fits well among the conventional supercon¬ 
ductors and sits close to Pb and Nb 3 Sn which are other 
highly optimized superconductors. Using this figure one 
can also translate to other figures for strong coupling 
corrections to BCS ratios as found in Refs—*— We have 
done the numerical calculation and find for the anhar¬ 
monic spectrum the specific heat ratio AC(T c )/"/T c = 
2.53 and the magnetic field ratio 7 T 2 /H 2 (0) = 0.136, 
which are very close to the values found for Pb, as ex¬ 
pected from Fig. IS2l 

We conclude from this that sulfur hydride appears 
to be a very highly optimized electron-phonon super¬ 
conductor which uses the high phonon frequency of the 
hydrogen and large electron-ion interaction to obtain 
a high T c . Nonetheless, it remains a question as to 
how to demonstrate the superconductivity and verify the 
electron-phonon interaction. Given that the sample is 
under pressure, many typical experimental probes may 
not be used. Tunneling was the classic experiment in 
conventional superconductors but here the energy scales 
are much higher than those that are usually accessed 
with tunneling. Consequently, an experiment which is 
non-contact but can probe high energy scales is needed. 
Here one is further helped by the fact that the mate¬ 
rial will be in the “clean limit”, where 2 Ao ^ 1 /r with 
1/r being the impurity scattering rate. 21 - An excellent 
probe in this case, then, is optical conductivity measure¬ 
ments as they may be done on a sample in a high pres¬ 
sure cell. As the cuprate superconductors have already 
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FIG. 2. (Color online) The gap ratio 2Ao/fcsT’c vs T c /uJi n . 
The solid black points are the data for several conventional 
superconductors taken from Table IV of Ref.—. Some have 
been labeled. The dashed line is an analytical formula devel¬ 
oped for conventional superconductors from strong coupling 
Eliashberg theory— The gold star is the value calculated for 
the anharmonic spectrum shown in Fig. ISlI and the red stars 
are the values for the harmonic spectra^ with T c = 190 K. 
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illustrated, the clean limit allows one to see boson struc¬ 
ture in optics.— This is illustrated in Fig. [S3] where we 
show our calculation for the optical conductivity for the 
anharmonic spectrum of H 3 S. [A comparison with that 
from the other a 2 F(u) spectra is given— in Fig. S2]. 
Shown are the superconducting and normal state curves 
for T = 0.1 T c . In the latter case, A is set to zero in the 
calculations! and this curve is shown at the same tem¬ 
perature to illustrate differences between the two states 
independent of temperature. [For the normal state at 
T c , which is not that different, see Fig. S3 of the Sup¬ 
plemental Material—] While the normal state shows a 
Drude absorption followed by a bump associated with 
the main source of inelastic scattering from the 200 meV 
peak in the electron-phonon spectrum (see Fig. ISTI). the 
superconducting state is much more dramatic. Here one 
sees a sharp impurity-dominated peak at the gap edge 
of 2Ao = 73 meV and then a series of absorption peaks. 
This structured absorption can be traced to the a 2 F{uj) 
spectrum which we show in Fig. [S3] shifted along the fre- 
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FIG. 3. (Color online) The real part of the optical conduc¬ 
tivity (scaled up by a factor of 10 4 ) vs oo at temperature 
T = 0.127,. Shown are the normal and superconducting state 
curves (red and blue, respectively). uj p is the plasma fre¬ 
quency. The a 2 F(uo) spectrum from Fig. IS1I is shown shifted 
to the right along the frequency axis by 2Ao. 

quency axis by 2A 0 . A full analysis and inversion— is 
needed to pull out the a 2 F(co) from the data, but it is 
clear that there are signatures of the energy gap, phonon 
structures and a shift in the main absorption hump going 
from the normal state to the superconducting state. By 
referring^ to Fig. S3 for temperature-dependent spectra, 
it can be seen that such structures and overall response 
survive up to very high temperature and are quite ro¬ 
bust. As a final note, for these calculations we have used 
a transport impurity scattering rate of 20 meV which is 
reduced by a 1 + A factor— to make the halfwidth on the 
Drude conductivity to be about 7.5 meV in this case. It 
should be remarked that the behavior seen here due to 
the clean limit is quite different from what is found for 
conventional superconductors, such as Pb^ which are in 
the dirty limit. The high energy scale of the supercon¬ 
ducting gap contributes to sulfur hydride being in the 
clean limit and showing strong Holstein structure. 

In summary, using three different electron-phonon 
spectral densities now available in the literature which 
are based on density functional band structure computa¬ 
tions for H 3 S under 200 GPa of pressure, we have verified 
that all give spectra that are near optimum for T c . More 
specifically all have a characteristic phonon energy wi n of 
Eliashberg theory which falls near the maximum of its 


own functional derivative of T c with respect to a 2 F(oj). 
Given that the functional derivative 5T c /5a 2 F(uj) dis¬ 
plays a very broad maximum, a significant distribution 
of phonon energies about wi n does not deplete T c much. 
For example, in the case of the a 2 F(co) provided by Errea 
et al> 3 including anharmonicity, the value of T c obtained 
falls only 8% below its theoretical absolute maximum for 
any distribution of phonon energies keeping the area un¬ 
der the spectral density fixed. 

In Eliashberg theory, it was found that a critical pa¬ 
rameter which measures coupling strength is the dimen¬ 
sionless ratio T c /u}\ n . For small values of this ratio BCS 
theory is recovered while strong coupling corrections to 
the universal dimensionless ratio of BCS increase with 
increasing value of T c /ui\ n . For H 3 S, one finds a value 
near 0.13 which is close to that found in Pb, Nb 3 Sn and 
others. It corresponds to the case where uq n is close to 
the maximum in the functional derivative T c /uj\ n ss 7 
and hence accounts for the near optimization of the spec¬ 
trum. This also implies that the thermodynamic prop¬ 
erties of H 3 S will be close to those of Pb and Nb 3 Sn. 
This does not imply, however, that other properties of 
their superconducting state will be similar. As an ex¬ 
ample, optical properties in H 3 S will be very different 
from those in other conventional superconductors which 
are in the dirty (impurity-dominated) limit— Instead, 
because of the very large energy scale associated with its 
phonons and its superconducting gap value, H 3 S is in the 
clean (phonon-dominated) regime. To be more precise, 
the absorptive part of the AC longitudinal conductivity 
Recr a . a ;(n) as a function of photon energy f 1 will show two 
distinct and separate regimes. At the energy of twice the 
gap, the regime remains impurity-dominated and there 
is a very sharp rise in absorption at S2 = 2A. but this is 
followed by dominant phonon-assisted absorption struc¬ 
tures extending to many times 2A which is usually re¬ 
ferred to as the clean limit regime because residual im¬ 
purity scattering plays no significant role. This region 
provides a detailed picture of the pairing glue. Conse¬ 
quently, optics can be used to provide information not 
only on the existence of superconductivity but also of 
the driving mechanism. 
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Shown in Fig.ISlIare three of the a 2 F(oj) spectra avail¬ 
able in the literature. The upper spectrum is a harmonic 
spectrum from Flores-Livas et at [Ref. 4] and the lower 
two are from Errea et al. [Ref. 3]. There is also a harmonic 
spectrum from Duan et at [Ref. 2] which is not shown 
here. These spectra were digitized from the original pa¬ 
pers and placed on a 2 meV grid which give character¬ 
istics that are in qualitative agreement with the original 
works if not always in perfect quantitative agreement. 
Such details do not impact our results and conclusions 
presented here. It is clear that all of the spectra are 
qualitatively similar. 



in the literature. The spectra have been offset for clarity. 

For our calculations, we have used the standard nu¬ 
merical solutions of the imaginary axis s-wave Eliashberg 
equations in order to solve for T c and other thermody¬ 
namic quantities. By an analytic continuation to the real 
axis, we evaluated the energy gap and the optical conduc¬ 
tivity. A description of these techniques and formulas can 
be found in the review by Carbotte [Ref. 9]. The optical 


conductivity calculation is done following the procedure 
outlined by Akis et al [Ref. 21]. 

In the following table we show the characteristic pa¬ 
rameters of our digitized spectra. 

In Fig. [S2] we show the conductivity spectra for each 


TABLE SI. Table of parameters for various digitized spectra 


a 2 F( u>) 

Errea 

anharmonic 

Errea 

harmonic 

Flores-Livas 

harmonic 

X 

1.67 

2.55 

2.47 

win (meV) 

122.1 

102.6 

106.1 

Wmax (meV) 

254 

218 

224 

A (meV) 

118.5 

146.6 

148.7 

Tc (K) 

190 

190 (249 )£ 

190 (255) 

h*(N = 6) 

0.18 

0.38 (0.16) 

0.40 (0.16) 

T C /UJ In 

0.134 

0.16 (0.209) 

0.154 (0.207) 

Win /T c 

7.46 

6.27 (4.79) 

6.48 (4.83) 

T c /A 

0.138 

0.112 (0.146) 

0.110 (0.148) 

Ao 

36.5 

39.9 

39.1 

2Ao/fcsT c 

4.46 

4.87 

4.78 


a Entries in brackets correspond to taking fi* = 0.16. 


of the a 2 F(co) shown in Fig. IS 11 To examine the effects 
of temperature, we show the conductivity for the anhar- 
monic spectrum at several temperatures in Fig. IS3I 



FIG. S2. (Color online) A comparison of the real part of 
the conductivity in the superconducting state for the three 
different spectra shown in Fig. IS 11 Here, T = 0.1T C . 
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FIG. S3. (Color online) Variation of the real part of the con¬ 
ductivity in the superconducting state for the anharmonic 
spectrum. Curves are shown for different temperatures as in¬ 
dicated. The red curve is for T = T c which is the normal 
state. 




